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Abstract 

For a piecewise linear version of the periodic map with anomalous diffusion, the 
evolution of statistical averages of a class of observables with respect to piece- 
wise constant initial densities is investigated and generalized eigenfunctions of the 
Frobenius-Perron operator are explicitly derived. The evolution of the averages is 
controlled by real eigenvalues as well as continuous spectra terminating at the unit 
circle. Appropriate scaling limits are shown to give a normal diffusion if the reduced 
map is in the stationary regime with normal fluctuations, a Levy flight if the re- 
duced map is in the stationary regime with Levy-type fluctuations and a transport 
of ballistic type if the reduced map is in the non-stationary regime. 
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1 INTRODUCTION 



The Frobenius-Perron (FP) operator is a powerful tool to investigate statisti- 
cal properties of dynamical systems [1,2, 3]. For strongly chaotic maps where all 
trajectories are exponentially unstable such as expanding and hyperbolic sys- 
tems, correlation decays exponentially and the decay rates, known as Pollicott- 
Ruelle resonances [4, 5], are given as the logarithms of the zeros of the dynamical 
zeta function. Note that the zeta function is essentially the Fredholm determi- 
nant of the FP operator and its zeros are eigenvalues of the FP operator in a 
generalized sense. Similar characterization was observed for periodic strongly 
chaotic maps with deterministic diffusion (see [2,3,6,7] and references therein): 
There, because of the periodicity of the system, the FP operator is decomposed 
into Fourier components and each component admits wavenumber-dependent 
eigenvalues including hydrodynamic modes of diffusion. Moreover, the corre- 
sponding eigenf unctions have fractality, which is closely related to the under- 
lying phase space structure and the transport coefficients [8]. 

The method of the dynamical zeta function is extended to certain non-hyperbolic 
maps with marginal fixed points[9,10]. For those maps, the power-law decay 
of correlation and/or the anormalous growth of the mean square displacement 
are obtained from the zeta function. The main difference in the zeta func- 
tion between the non-hyperbolic and hyperbolic maps is the existence of a 
cut terminating at the unit circle in the former case. In other words, the FP 
operator of the non-hyperbolic map has a continuous spectrum terminating at 
the unit circle. However, its spectral structures such as the eigenfunctions as- 
sociated with the continuous spectrum are not fully understood. Recently, as 
a first step towards the full spectral characterization of intermittent chaos in 
non-hyperbolic maps, we have obtained generalized eigenfunctions of the FP 
operatorfll] for a piecewise linear version of the Pomeau-Manneville map [12] 
introduced by Artuso[10]. We continue the study of Ref.fll] for non-hyperbolic 
maps with anomalous diffusion and, in this paper, we investigate generalized 
eigenvalues and the associated eigenfunctions of the FP operator for a piece- 
wise linear version of the map studied by Geisel, Nierwetberg and Zacherl[13]. 
The mean square displacement of this map spreads faster than the normal 
diffusion, or it exhibits super-diffusion. Maps with sub-diffusion such as that 
studied by Geisel and Thomae[14] will be discussed elsewhere. 

Note that anomalous transports such as the Levy diffusion and Levy flight are 
known to appear in a variety of natural systems and have been extensively 
studied (see [15] and refereces therein), and still new maps with anomalous 
diffusion are proposed[16,17]. Then, it is worthwhile to develop a dynamical 
systems approach because it would provide a microscopic understanding of 
the anomalous transports. This is one of the objectives of this work. 
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The paper is organized as follows. In Sec. 2, we describe a model and classes 
of observables and initial densities. Evolution of statistical averages of ob- 
servables is investigated in Sec. 3 and generalized eigenf unctions of the FP 
operator are derived. The behavior of the wavenumber- dependent eigenvalues 
responsible for the anomalous diffusion is investigated in Sec. 4. As one of the 
applications, 'macroscopic' transports are considered in appropriate scaling 
limits. The normal diffusion is derived if the reduced map is in the stationary 
regime with normal fluctuations and Levy flights are obtained if the reduced 
map is in the stationary regime with Levy-type fluctuations. And if the re- 
duced map is in the non-stationar regime, a new transport equation is shown 
to be obtained. Conclusions are drawn in Sec. 5. 



2 Piecewise Linear Geisel-Nierwetberg-Zacherl Map 



One of the simplest maps with super-diffusion is that introduced by Geisel, 
Nierwetberg and Zacherl[13]. Here we consider its piecewise linear version <fi. 
The map is defined on the whole real line (Fig.l (a)): : R — > R as 



(f>(x) 



<f>(x - [x]) + [x] - 1 , if x - [x] < 1/2 
4>(x - [x]) + [x] + 1 . if x - [x] > 1/2 



The symbol [x] stands for the smallest integer which does not exceed x, and 
4> : [0, 1] — > [0, 1] is the reduced map (cf. Fig.l (b)) defined as follows: 



<t>{y) 



rik{y-ik) + i 



k-1 , 



(6 < 2/ < 6-0 



Vk(y - i + 6) + i - 6-1 , (i - 6-i < y < 1 - 6 



(2) 



where k — 1, 2, • • • and the variables 6 an d i]k are given by 



Vk = < 



2C(£) 

6-2 — 6-i _ [ _ 
6-i — 6 \^ 



' 6 "2 



I 1-2^ 



—J . (fc = 2,3,..-) 
(k=l) 



(3) 



(4) 



In the above, (3 > 1 is a parameter and £(s) = S^Li n ~ s is the Riemann 
zeta function. The map has marginal fixed points y = and y = 1 where 
the map is approximated, respectively, as <p(y) — y + Coy /3 ^ /3 " 1 - ) and <f)(y) — 
y — c (l — yy/(P~ 1 } with a positive constant c . When (3 > 2, the map is 
in the stationary regime and, when /3 < 2, it is in the non-stationary regime. 
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2 
1 



■1 

-1 



x 




(a) 

Fig. 1. (a) The piecewise linear Geisel-Nierwetberg-Zacherl map <j) and (b) its re- 
duced map 4> (/? = 2). 

Moreover, fluctuations in the stationary regime are normal when (3 > 3 and of 
Levy-type when 3 > j3 > 2 [18]. Those fluctuations lead to anomalous diffusion. 
Note that a similar piecewise linear map was studied by Artuso[10]. 

The FP operator P and its adjoint P* are defined as 



Pp(x) = j dyb (x - p(y) 

— oo 

P*A(x)=A(J)(x)) . 



(5) 
(6) 



Thanks to the periodicity of the system, it is convenient to introduce the 
Fourier representation with respect to the integer part of the coordinate [x] : 



A(q,y)= J2 e^A{n + y), p(q,y) = £ e in «p(n + y) , 



(7) 



where y G [0, 1). Then, the average value of an observable A at time t with 
respect to an initial density p is given by 



+7T 1 

(A,P t p)^(P* t A,p) = f£f dyMq,yyp(q,y) 

-IT 



(8) 



where (A, p) is the inner product: 



11 dq 1 

(A./)) = j dxA(x)p(x) = J ^ J dyA(q,y)*p(q,y) 

-IT 



(9) 
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and A t (q, y)* is the Fourier transform of P U A obeying the recursion relation: 



In order to extract information about generalized eigenfunctions of the FP 
operator, we consider the time evolution of the averages of observables with 
respect to a class of initial densities. 

The Fourier transform of observables A(q, y) are assumed to be smooth at 
the origin and at y — 1 in the sense that they behave as A(q,y) = A(q, 0) + 
a i(q)y + O (y 3 /^- 1 )^ near y = and as A(q,y) = A(q, 1) + a 2 (q)(y — 1) + 
near y — 1, where ai and a 2 correspond, respectively, to 

d y A(q, y)\ y =o, and d y A(q,y)\ y=1 . More precisely, is a bounded function 

such that there exist functions ai(q) and a 2 (q) of g and inequalities 



\A(q,y)-A(q,0)-a 1 (q)y\<K\y\^^ (11) 
y) - A(q, 1) - a 2 ( 9 )(y - 1)| < K'\l - yf^ (12) 

hold for < y < 1 and — 7r < q < tt, where K and K 1 are positive constants. 
Let Xo be a set of such functions, then it is a Banach space with respect to 
the norm: 



ii /. I, v-* in \A{q,y) - A(q,0) - a l (q)y\ 

UWo = E 2 Halloo + Woo + jup wmzT) 

Q11 „ \A(q,y)-A(q,l)-a 2 (q)(y-l)\ 
q ,o<y<i |1 -y\ p '^ L > 

where || - || oo = su Pg,x I ' I- The space X is invariant with respect to the 
adjoint of the FP operator P*: namely, if A e X , then P*A e X - Note that 
since X contains a space C 2 [—ir, it] x [0, 1] of twice continuously differentiable 
functions and C 2 [— it, n] x [0, 1] is dense in the Hilbert space L 2 [— it, n] x [0, 1] 
of square integrable functions, X Q is dense in L 2 [—-n,-n] x [0, 1]. 

It was shown [4,5] that the spectrum of the FP operator of expanding maps 
is controlled by the smoothness of densities. In this proof, the positivity of 
the Lyapunov exponents plays an essential role. However, the map <fi does 
not enjoy such a property and the general situation is not known. So, we 
restrict ourselves to a narrow class of piecewise constant initial densities in 
order to derive explicit expressions of the generalized eigenfunctions of the FP 
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operator. More precisely, we consider a set Xd of initial densities p given by 

p{q,y) ~\pt(<i), <y< 1-6) ( j 



where A: = 1, 2, • • - and 

+ 0O / h, \ P +00 

^) = ^pHq)[ 1 — i j with EII« t M'<+ 00 ( 15 ) 



for some constant 6 > 1. The space X D is a Banach space with respect to the 
norm 

D = Y^{\\pt Woo + \\pt \\oo}e l , (16) 

1=0 



and is invariant under the action of the FP operator P: namely, if p G X D , 
then Pp E X D . Note that, in contrast to the space X Q of observables, the 
space X D of densities is not dense in the Hilbert space L 2 [—n,n] x [0, 1] of 
square integrable functions. 



3 Evolution of statistical averages of observables 

3. 1 Matrix elements of the resolvent 



In this section, we investigate the evolution of the average (A, P l p) of an 
observable A with respect to a piecewise constant initial density p. We begin 
with the Neumann series of the matrix element of the resolvent: 



1 + °° 1 

(A--^p) = E^tt(ApV) (17) 



where 1 is the identity operator. In terms of 



3 k (t,q)= J dyA t (q, y y , and B+(t,q)= J dyA t (q,y)* , (18) 



I'" 

each term of the series can be rewritten as 
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(AP'P) = I ^-T,{pt(Q)BUt,q)+P^q)Bj:(t, q )} , (19) 

_,r Z7F k=0 

and, hence, 

1 + r Ha +oc 

(A,—~P) = J ^T,{p^9)Bt(q,z)+p^q)B^q,z)} , (20) 



where B^(z,q) = E t = (*, q)/z t+1 . The right-hand side of (20) is calcu- 
lated as follows: For example, when k > 2, eq.(2) gives recursion relations of 



£? fc (t + l,g) = J dye~ i «A t (q,ri k (y-Zk)+Zk-iY 

-ig f ^ -iq 

= / dyA t (q,yy = B^q), (21) 

Vk •> Vk 
Sfc-1 

e iq 

B+(t + l,q) = —B+_i(t,q) , (22) 
Vk 

and, hence, we obtain 

1 e ±iq - 

B±(q, z) = -B±(q, 0) + —B}_i(q, z) (23) 

Z ZT] k 



Similarly, one has 

1 p ±i( i +°° 

Bf(q, z) = -Bt(q, 0) + — £ *) (24) 



The recursion relations (23) and (24) give 



2^Pk(q) B k{q,z) = = + 

fc=i z 



+ wlw* (9 ' 2) + Jcw* (? '^ (25) 



where S ± (g,z), $ ± (g, z), \l/ ± (g,z), f2(g, z) and Z(z) are defined by: 
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1=2 Z I k=l k 

k / 7 \ Z 3 D± 



(26) 



CO 



fc =l ^ 



V(q,z) = l 2 



(29) 

= E ^p- (3°) 



fc=l 



Infinite series in (26), (27), (28) and (30) are absolutely convergent for \z\ > 1 
and, thus, are analytic there. In addition, because of (17) and the analyticity 
of the resolvent in \z\ > 1, one finds 

|z|=r — T z|=r 



where R(q,z) = J2T=i{pt (q)B£(q, z) + p k (q)B k (q, z)} and the integration 
path is a counter-clockwise circle centered at z = with radius r(> 1). 



3.2 Analytical properties of individual functions 



Here we derive analytical continuations of $ ± (g, z), \l/ ± (g, z) and Z(;z) into 
the unit disk \z\ < 1. With the aid of the formula r(/5)/^ = / °° dss^e"**, 
the series expression of $ ± can be rewritten as 



••—S^/---' 5(5=)' 

^ e ±i( iV((3) J ze^v-e- 8 1 ; 



o 

The above calculations are justified because the summation converges uni- 
formly in s provided \z\ > 1. 

We observe that, as a result of an inequality for bounded observables A: 
l p \Bf(q,0)\ < -6)ll^l|oo = ||A|| 00 /(2C(/3)), the estimate 
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E 



lf>B±(q,0) 



ds 



s P-l e -U 



ze ^q 



< 



\A\ 



d(ze^i, [0,1]) 



(33) 



holds where d(£, [0, 1]) is the distance between £ and the real interval [0, 1]. It 
is obvious that each term of ^ in the right-hand side of (32) is analytic except 
a cut on an interval {z = te ±tq : < t < 1}. Therefore, (32) defines a function 
which is analytic except the cut and, then, is an analytical continuation of $ d= 
from the outside of the unit disk to the whole complex plane. 

Similarly, the following analytical continuations are obtained: 



oc 



u 



P-l p -(l+l)u 



ze Tn — e - 



ds 



(34) 
(35) 



both of which are analytic in the unit disk \ z\ < 1 except the cuts. The function 
\l/ ± (g, z) has a cut on an interval {z = te ±iq : < t < 1} and Z(z) has it on a 
real interval [0, 1]. In addition, the first two terms of (25) have the following 
analytical continuation: 



Z " 1=1 



The series of the right-hand side are absolutely convergent for \z\ < 1 and, 
thus, analytic there. Note that (36) is not a simple rewrite since the series in 
(26) does not converge for \z\ < 1. 

As a result, the integrand of (31) is meromorphic in the unit circle except 
the cuts. Its poles, if they exist, are given as zeros of Q(q,z). When q = 0, 
Q(0, z) = is equivalent to Z(z)/ z±(((3) = 0. As easily seen, Z(z)/z — ((j3) = 
has a unique zero at z — 1 and an argument similar to that given in Ref . [1 1] 
shows that Z(z)/z + ((f3) = has a unique real zero on (—1,0). Or Q(0,z) has 
two real zeros; one is in the interval [—1,0] and the other at 1. When q ^ 0, 
the behavior of Q(q, z) is investigated numerically with the aid of Newton's 
method and we found that f2(g, z) may have, at most, two simple real zeros; 
one in the interval [—1, 0] and the other in [0, 1], which we denote, respectively, 
as \j(q) and Aj(g). Their properties will be discussed in Sec. 4 again. 
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3.3 Decomposition of averages of observables 

Now we go back to the contour integral in (31). Note that, because of the 
g-integration, it is sufficient to consider g ^ terms. Then, the arguments 
of the previous subsection show that, for each q, the integrand is analytic 
in \z\ < 1 except the simple poles at z = Aj~(g) and cuts on the segments 
{z = te ±iq : < t < 1}. Therefore, the contour integral is evaluated by 
deforming the contour from the circle \z\ — r to the counter-clockwise curve 
C(q) shown in Fig. 2 and one gets 

f ^- jR(q, *) = £ Resxr/Riq, z) + f ^ jR(q, z) (37) 

\z\=r a=± C(q) 

where X^(q) are abbreviated as A^ and Res A ± stands for the residue at z = A^ . 

d 

Firstly, we consider the contribution from the curve C(q) in case of q > 0, 
which consists of two counter-clockwise arcs ^or = {z '■ \z\ = e & |argz| < \q\} 
and 7ol = {z : \z\ = e & \ir — arg^| < n — \q\}, two counter-clockwise circles 
7± = {z : \z — e ±ig \ = e}, the segment C e +l = {z = (t — i0)e tq : e < t < 1 — e} 
(C e +U = {z = (t + i0)e iq : e < t < 1 - e}) just below (above) the up- 
per cut and the segment C e _i = {z = (t — i0)e~ ig : e < t < 1 — e} 
(Ct u = {z = (t + i0)e~ iq : e < t < 1 — e}) just below (above) the lower cut: 

C(q) = Ior U C £ +l U 7+ U C\ u U 7ol U C e _ x U 7 _ U CI U . 

Here, e stands for a small positive number, C^_ u and Ct u point to the origin, 
and C e +l and Ci^ point to the infinity. 




Fig. 2. Integration contours (solid curves) for individual q. Two dots and dotted 
lines indicate, respectively, the simple poles A . (q) and the cuts of the integrand, 
(a) The original circle \z\ = r(> 1) and (b) the deformed contour consisting of two 
circles and the curve C(q) surrounding two poles and the cuts, respectively. The 
directions are counter-clockwise and the unit circle is shown as a reference. 
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Since the conditions A e X and p G lead to 

lim\z\ a \R(q, z)\ = (0 < 3a < 1) , 

the contributions from 7 _r and 7 l vanish in the limit of e — > 0. Similarly, the 
contributions from 7± vanish in the same limit since 

lim \z-e ±iq \ a '\R(q,z)\=0 (min(2 - /3, 0) < a' < 1) . 
On the other hand, we have 



(38) 



where = = te* 9 : < t < 1} is the upper cut and z+m = z + (— )ie* 9 0. 
Eqs.(25), (32) (34) and (35) lead to 



R(q,z+)-R(q,z+) 



2m 



= N + (q,z) 



A^ + (q,z) + 



ze*C(P) 



xze 



-iq 



(39) 



where 



N + (q,z) = 
A^ + (q,z) = 
A$ + (q,z) = 



ln(ze~ 



iq 



1/3-1 



r(/3) [fi(g, <)^ + (g, *, + ) - ^ + )^ + (g, z' 



2mze- i( i [- ln^e"^)]^ 1 

T{0) [Q(q, 4)<t> + (g, - n((z, ^ + )$ + (g, 4) 



27rize- < « [- ln(ze"^)] 



(40) 
(41) 
(42) 



The contributions from the cut with g < as well as from the other cut 
C~ = {z = te~ iq : < t < 1} are evaluated in a similar way and finally one 
obtains 



lim 



/ 2% ztR ( q ' Z "> = J dz (4*±(?,*)) zt (F±{q,z),p) (43) 



The functional (A, F±(q, zfj is defined by 
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1 



(A,F ± (q,z)) = J dyA{q,yyg%{y) , (44) 
o 

+E'Vfa) ^ )2r(/ ; ) /^p^ ¥ ^. (45) 

xTto) = rf(i-») = {J J^iJ'- 1 ' • («) 

! - ^ L s P^£^ , (47) 



where P stands for Cauchy's principal value. Note that ze Tiq is real if z e C^. 
Similarly, (^F±(q,z),pj is given by 

pf(g)Z(ze ±i,? ) 7 J s /?-i e -0+i) S 



n^g^K^ + Scfe /^^Fr} ' (48) 



92 r(/3)n(q,z + ie ±iq O)n(q, z - ie ±iq O) ' 1 ' 

The pole contributions are calculated as above and one has 



E Res A ^i?(g, z) = £ A5(9)' lim (* - X d (q)) i2(g, *) 
<j=± a=± z —* d yi) 

= J2{AFZ( q )) K( q y ( F d(q),p) , (50) 

(7=± 

where the functionals are given by 



i 

(A, F°(q)) = J dyA(q, y)*g±(q, y) , (51) 
o 

OO 



+ fe^X^/" S A^-e-' (52) 
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e"'V, + (<j) 7. s"-'e-« +1 )' 



(^M,p)=g fi , ( ; A ^ ( ; )r(/ , ) / rf ^ 

oc 



± e~ ic i - e~ s 



In the above, A^(g) is abbreviated as A J and Cl'(q, z) stands for the derivative 
of Q(q, z) with respect to z. 

In summary, we have derived 
dq 



(A,P t p) = J g E(A^)) A5( 9 )* (^(g),p) 

— 7T O"— i 

+ E (A^M) ** (^(g,^),p) (54) 



= lt C Q 



3.4 Generalized eigenf unctions of the Frobenius- Perron operator 



We remark that, for any bounded linear functional F over Xo, (P*A, F) (A e 
X G ) is well defined due to the invariance of X Q under P* and that the operator 
P can be extended via (A, PF) = (PM., F) to the space of all bounded linear 
functional over X (the dual space of X )- The linear functionals F^ and 
F±(q, z) are bounded and satisfy the following relations for any A e X : 



(A P F d(q)) = (P*A, Fj(q)) = \a(q)(A, F^(q)) , (55) 
(A,PF±(q,z)) = (P*A,F±(q,z))=z(A,F ± (q,z)) . (z e Cf) (56) 

Namely, they are eigenfunctions of P in the generalized sense [19]. 

Similarly, the operator P* can be extended to the dual space of X D and the 
functionals F^, Fi(q,z), F u (q,z) are its generalized eigenfunctions: 



(P*F±(q),p) = (F±(q),Pp) = )${q)(F?(q),p) , (57) 
(P*F ± (q,z),p) = (F ± (q,z),Pp) = z(F ± (q,z),p) , (z e Cf) (58) 

Sometimes, eigenfunctions of P and P* are referred to, respectively, as right 
and left eigenfunctions of P. 

It is then interesting to reinterpret (54) from this point of view. One has 
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(Ap) = /^E(A^(?))(^(?),p) 

+ /^E J dz(A,F a (q,z))(F a (q,z),p) (59) 

-7T Q = ± C« 

+ /^ E (A^M) « (60) 



The first relation (59) corresponds to the completeness of the set of left eigen- 
functions {F^(q), F±(q, z)} and that of right ones {F^(q), F±(q, z)}. The sec- 
ond relation (60) can be considered as a spectral decomposition of P in terms 
of those eigenfunctions, which correspond, for each q, to two isolated eigenval- 
ues \a(q) and continuous spectra on C^. These relations are regarded as an 
extension to the FP operator of the generalized spectral decomposition in the 
sense of Ref. [19], which was originally formulated for self-adjoint and unitary 
operators. Such extended spectral decompositions are widely used to study 
statistical properties of dynamical systems (see [2,7] and references therein). 



4 Transports 



For systems exhibiting deterministic diffusion with diffusion coefficient D [2,3,6,7], 
the FP operator has an eigenvalue z — z q such that z q ~ 1 — Dq 2 for small 
wavenumber \q\ and it fully characterizes the diffusion in a scaling limit. A 
similar behavior is expected for the map <j). We then investigate an eigenvalue 
A J (q) for \q\ <C 1 and appropriate scaling limits. 



4.1 Case of (3 > 3 



When \q\ <C 1, |1 — \^e ±iq \ <C 1 and it is sufficient to retain the largest terms 
of £l(q, z) with respect to 1 1 — ze ±iq \. For f3 > 3, they are found to be 



Sl(q, z) ~ ^^ {zcosq- 1) - (^f) (1 + 2^cosg) 
C(/3-l) + C(/3-2) /i o 2 



{l — 2zcosq + z cos2g) . (61) 
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Fig. 3. The ^-dependence of the eigenvalues of the FP operator. The upper branch 
is At(g) and the lower branch is AT(g). (a) f3 = 4, (b) /3 = 2.5 and (c) /? = 1.7. 

Then, the equation fi(g, A^) = has a solution 



A d — 



1 

1 - - 

2 



CQ9-2) C(/3-l) 
CC9-1) C(/3) 



s 2 , 



(62) 



which implies that the map exhibits ordinary diffusion. In Fig. 3 (a), the q- 
dependence of numerically calculated X d (q) for (3 = 4 is shown. The q 2 - 
behavior of A^(g) near q = is clearly seen. 

The appearance of ordinary diffusion can be seen more explicitly by consid- 
ering a scaling limit, where one sets n = L^X, t = L T T and the limit of 
— > +oo is taken while assuming a functional relation between and L T . 
Accordingly, the observable A(L^X) is assumed to be a well-behaved function 
of X. Then, one has 



1 - 1 oo 

— A(Q/L 5 ) = — e inQ / L tA(n)~A^\Q) 



(63) 



where A^ cg) (Q) = J^ 00 dXe l( ^ x A(L^X) is the Fourier transform of the coarse- 
grained observable. Similarly, the Fourier transform of the probability distri- 
bution is assumed to be constant on each unit interval and is a slowly varying 
function of the site coordinate [x] . Then, in terms of the Fourier transform of 
the coarse-grained density p( cg \ one has pf{Q/L^) = Si i0 p^ cg '(Q). Hence, 



+ Yl e m<?( ^ T+1) Ids (e, F a (q, se aiq )) s L ^ T (F a (q, se mq ), e)\ (64) 

where e is a function of n and y with the Fourier transform e(q,y) = 1. In 
order to have a non-trivial scaling limit, one must have L T oc L? and here we 
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take L T — L|. Then, one has 



-Dg Q T 



(65) 



L T =Li 



where Dp 



N _ 1 



C(/3-2) _ CQ9-1) 
C09-1) C(/3) 



is the diffusion coefficient. In this scaling limit, 



the contribution from \ d vanishes as sup^ \\ d (g)| < 1. This is also the case for 
the continuous spectrum as k(s, Q) = lim (e, F a (q, se aiq ))(F a (q, se aiq ), e)\ q_ 



is a bounded function of s: 



lim 



jds (e,F a (%,se aiq )) s L ^ T (P a se aiq ) , e) 
o 

i 

< sup|fc(s,Q)| lim [dss LrT = 0. 



(66) 



In short, the scaling limit of the average value (A) T = limi^oo (A,P LtT p) of 



A at time T is given by 
dQ 



(A) T = J ^A^iQY e~ D ^ 2T p^(Q), 



(67) 



which is nothing but the solution of the diffusion equation. 



4.2 Case of 3 > (3 > 2 



When (3 = 3, the largest terms of £l(q, z) with respect to j 1 — ze ±iq \ are 



n , v 2C(2) (l-ze^) 2 log(l-^) 

fi(^) - ^y(^cosg- 1) + ^ , (68) 



from which one obtains the following solution of £l(q, Aj~) = 0: 
4C(2) 

Similarly, when 2 < (3 < 3, the leading term of Vl(q, z) is 
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^-^(^-1)+^ f ^ y — 2j . (70) 
and the approximate solution of fi(g, A J) =0 is 

(71) 



7r \q 


0-1 


2C(/3 - i)r(/3) 


cosf 



In both cases, the system exhibits anomalous diffusion. The g-dependence of 
numerically calculated A^ for (3 = 2.5 is shown in Fig.3 (b) and a g-dependence 
different from that for the normal diffusion is seen in the figure. 

Now we consider a scaling limit for 2 < (3 < 3, where we set n = L^X, 
t = L T T and take the limit of — > +oo while keeping L T = L^~ l . Then, 
the contributions from the eigenvalue A^ and continuous spectrum vanish 
because of the same reason as in the previous subsection, and in terms of the 
coarse-grained observable A^\Q) and probability density p^ c% \Q) introduced 
previously, one has 

{A) T = lim {A,P L * T p) = J ^-A^\QT e' D ^- lT p™(Q) , (72) 



4.3 Case of 2 > 13 > I 

When (3 = 2, the leading term of £l(q, z) with respect to 1 1 — ze ±iq \ is 

n , , l-ze iq , l-ze iq l-ze~ iq , 1 - ze~ iq 

fife, log — - + -^- hg — ^. (73) 

By setting -ze~ iq = e i7T + re id and -ze iq = e in + re i{2n ~ e) (0 < 9 < 2tt), one 
finds the approximate solution of Q(q, Aj) = 0: 

A ^ 1 "OT' (l,l<<1) (74) 

On the other hand, when 1 < f3 < 2, the leading term of fi(g, z) is 
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nM * < (1 - ze ' r " + (1 - 2e "'* ) '" 1 } ' (75) 

and, by the same way as above, one finds that fl(q, Aj) = admits a solution 



q\ , (\q\ « 1) (76) 



when 3/2 < (3 < 2 and has no solution when (3 < 3/2. The g-dependence of 
numerically calculated A^ for (3 = 1.7 is shown in Fig. 3 (c), where |g|-linear 
behavior of A J is seen near q = 0. Moreover, for /3 = 1.1, Newton's method 
fails to find a solution of z) = and this is consistent with the absence 
of the solution. 

Finally, we consider a scaling limit for 1 < f3 < 2, where we set n = L^X, 
t = L T T and take the limit of — > +oo while keeping L T = L^. As be- 
fore, the contribution from the eigenvalue A^ vanishes. In this case, however, 
the continuous spectra do contribute and, after a tedious but straightforward 
calculation, one obtains 



8111 V2p=l)J 



sin ( 



l\) 



cot 



7T 



209 -1) 



OO 

(A) T = lim (A,P^ T p)= f ^-A^XQ)* pte\Q) (cosQT 

J 2,71 I 

L T = L^ — OO 

r rfz e-IQI rz (g - O^ 3 - 1 - (g + \ 

+ J 2tt z 2 + i (z-iy-i + iz + iy- 1 ) ' ( ^ 

where r = = t + i(l + e)|t : +oo -> 0} U = (1 + e)e if? |# : f -> ^}U 
{2; = t — «(1 + e) 1 1 : — > +00} is the integration contour with e an arbitrary 
positive constant, and the branch cut of (z =F i) 13 " 1 is {2 = t ± «|t : — > +00}. 



5 CONCLUSIONS 



We have studied, for a piecewise linear approximation of Geisel-Nierwetberg- 
Zacherl map, the evolution of statistical averages of bounded observables which 
are differentiable at the origin with respect to piecewise constant initial densi- 
ties which have certain smoothness at the origin. Because of the periodicity of 
the system, the average value is expressed by a superposition of Fourier com- 
ponents and the time evolution of each component is governed by two simple 
eigenvalues X^(q) and continuous spectra = {z = te ±iq : < t < 1} of 
the FP operator. The corresponding generalized eigenfunctions are explicitly 
constructed. 
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The main difference from the hyperbolic systems with deterministic diffusion 
is the existence of the continuous spectra starting from the unit circle for 
each wavenumber q. In spite of this fact, when /3 > 2, the transport in an 
appropriate scaling limit is entirely controlled by a simple eigenvalue X^(q) 
located near 1 and the average value of the coarse-grained description (A) T 
obeys 

(A) T = J ^AM(QY e' D ^ T p^\Q) , (78) 

— oo 



where \x = 2 and Dp = Dp for p > 3, and /i — P — 1 and Dp = Dp 
for 3 > (3 > 2. In other words, when the reduced map is in the stationary 
regime, transports in the scaling limits are governed by the eigenvalues. And 
the normal diffusion appears if the fluctuations are normal and the Levy flight 
appears if the fluctuations are of Levy-type. 



Here a remark is in order. As is well known, in periodic maps with marginal 
fixed points of Geisel-Nierwetberg-Zacherl type, the mean square displacement 
is well-defined and grows faster than the normal diffusion. Indeed, with the 
aid of (54), one recovers a known result on the mean square displacement: 



<(M - <M)t 



)t oc 



t 

tfogt 
t 2 



P>3 
P = 3 
3 > > 2 
2 > p > 1 



(79) 



where (• • -) t is the average at time t with respect to an initial distribution 
p{x) = 5[x],ra f° r some no- On the other hand, the mean square displacement is 
ill-defined for Levy flights[20]. However, there is no contradiction between the 
two conclusions because the ill-definedness of the mean square displacement 
comes from the scaling limit. 

The situation changed for 2 > f3 > 1 or when the reduced map is in the 
non-stationary regime. Then the 'macroscopic' transport is described by (77), 
the first term of which represents ballistic transport, while the second term 
corresponds to a more complicated process. 

Note that the space Xp of initial densities is not dense in L 2 [0, 1] and we will 
discuss elsewhere a complementary description in terms of initial densities 
which form a dense subset of L 2 [0, 1]. 
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